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ABSTRACT 

Prominent Fc Ka line emission is detected in the XMM-Newton spectrum of the 
ultraluminous infrared galaxy Arp220. The centroid of the line is found at an energy 
of 6.7 keV and the equivalent width of the line is EW ~ 1.9 keV (at 3.5c significance). 
A few other spectral features are found at various degrees of significance in the lower 
energy range on a hard 2.5-10 keV continuum (r ~ 1). The large EW of the Fe K line 
poses a problem with interpreting the hard X-ray emission as integrated X-ray binary 
emission. A thermal emission spectrum with a temperature of kT ~ 7 keV modified 
by absorption of Ah — 3 x 10 22 cm -2 , can describe the 2.5-10 keV continuum shape 
and the Fe K emission. A hot bubble that is shocked internally in a starburst region 
would have a similar temperature and gives a good explanation for the observed X-ray 
properties with a high star formation rate. An ensemble of radio supernovae in a dense 
environment, as suggested from VLBI imaging, could be another possibility, if such 
powerful supernovae are produced continuously at a high rate. However, the apparent 
lack of emission from X-ray binaries is incompatible with the high supernova rate 
(~ 2 SNe yr" 1 ) required by both interpretations. Highly photoionized, low-density 
gas illuminated by a hidden Compton-thick AGN is a possible alternative for the hard 
X-ray emission, which can be tested by examining whether radiative recombination 
continua from highly ionized Ca and Fe are present in better quality data from a 
forthcoming observation. 
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1 INTRODUCTION 

Arp220 is the nearest ultraluminous infrared galaxy (ULIG) , 
at 2 = 0.018 (corresponding to d = 74 Mpc for Ho = 
75km s" 1 Mpc" 1 ), with > 95 per cent of its total bolo- 
metric luminosity emitted at infrared/submillimetre wave- 
lengths [L ir (8-1000^m) ~ 1.3 X 10 12 L©] (Soifer et al. 1984; 
Sanders et al 1988, 2003). It is an advanced merger system 
involving two relatively large spiral disks, as evidenced by 
the detection in the optical of two large, faint, crossed tidal 
tails (Joseph & Wright 1985; Sanders et al 1988), with two 
nuclei separated by ~ 1 arcsec (~ 300 pc), as determined 
from images in the radio (e.g. Becklin & Wynn- Williams 
1987; Norris 1988; Baan & Haschick 1995) and near-infrared 
(e.g., Graham et al 1990; Scoville et al 2000) bands. High 
resolution CO observations imply an extreme nuclear molec- 
ular gas concentration (~ 3 x 1O 9 M0 at radii < 300pc cor- 
responding to a mean H 2 density of 10 4 -10 5 cm -3 ) (e.g., 
Scoville et al 1986; Solomon, Radford, Downes 1990; Scov- 



ille et al 1998; Sakamoto et al 1999). This is probably a 
consequence of disk gas being funnelled into the central re- 
gions of both galaxies during the merger (e.g., Barnes & 
Hernquist 1992). There is evidence for a galactic-scale out- 
flow, as Ha and soft X-ray nebulae with a size of ~ 30 kpc 
(Armus, Heckman & Miley 1990; Heckman et al 1996) sug- 
gest, which is thought to be driven by a super-starburst (but 
see Colina, Arribas & Clements 2004 for a merger shock in- 
terpretaion). The optical spectrum is LINER- like (Veilleux, 
Kim & Sanders 1999), most likely due to shock-heated inter- 
stellar gas by the galactic wind (e.g., Taniguchi et al 1999), 
while the double nucleus region is heavily obscured by dust. 

Arp220 is often used as a nearby template of lumi- 
nous star-forming galaxies at high redshift, thus verifying 
the true dominant power source in Arp220 continues to be 
of great importance. There is still no convincing direct evi- 
dence, from radio to hard X-ray wavelengths, for an active 
nucleus in Arp220. The general assumption has been that 
the bolometric luminosity of Arp220 is mostly powered by 
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a starburst, as perhaps best represented by a view based 
on ISO mid-infrared spectral characteristics (e.g. the PAH 
feature at 7.7 /im - Genzel et al 1998; Sturm et al 1996). 
However, the interpretation of the ISO data for ULIGs, and 
Arp220 in particular, has recently been revised by Spoon 
et al (2004), who took ice absorption into account and con- 
cluded that the PAH component is rather weak with only 
moderate obscuration. They inferred that a major power 
source, whether a super star cluster or an AGN, must be 
deeply enshrouded in dust, similar to an interpretation pre- 
viously proposed by Dudley & Wynn- Williams (1997). A few 
authors have suggested that a hidden QSO could be a major 
power source for Arp220 (e.g., Sanders et al 1988; Haas et 
al 2001). If an energetically significant AGN is present in 
Arp220, it must be Compton-thick (e.g., Rieke 1988). Hard 
X-rays constrain the lower bound of the absorbing column 
density to be Ah ~ 10 25 cm -2 (Iwasawa et al 2001). The 
lack of cold reflection characterised by a 6.4 keV Fe K line 
means that the covering factor of the obscuring matter has 
to be close to unity. 

Recent studies of nearby elliptical galaxies suggests 
that all large spheroids contain a supermassive black hole 
(SMBH: e.g., Richstone et al 1998). Given that the host 
galaxy of Arp220 already appears to have relaxed into an 
elliptical-like r 1/4 -law profile (e.g. Wright et al 1990), it 
seems reasonable to speculate that Arp220 may also contain 
a SMBH, perhaps even two, since the two observed nuclei 
have yet to merge. Two SMBH have recently been discov- 
ered in the luminous infrared galaxy NGC6240 (Komossa et 
al 2003), which is also an advanced merger system similar to 
Arp220. Given the detection of an enormous nuclear concen- 
tration of molecular gas in Arp220 there is plenty of material 
in the nuclear region to feed, and therefore build a SMBH, 
while also simultaneously suppressing the amount of direct 
emitted radiation from the nuclear source. This is not an 
implausible scenario, and a sensitive hard X-ray observation 
would have a better chance to catch the faintest sign of such 
a hidden AGN than observations at longer wavelengths, as 
sensitive hard X-ray observations of NGC4945 (Iwasawa et 
al 1993; Done et al 1996; Guainazzi et al 2001) and NGC6240 
(Iwasawa & Comastri 1998; Vignati et al 1999; Ikebe et al 
2000) have shown. 

We present here our analysis and interpretations of the 
public archive XMM-Newton data on Arp220, with which 
a strong Fe K line is detected for the first time. This will 
hopefully be followed by a longer observation that should be 
scheduled sometime during AO-3. 



2 OBSERVATIONS AND DATA REDUCTION 

Arp220 was observed with XMM-Newton on 2002, August 
11, and 2003, January 15. The two XMM-Newton observa- 
tions were carried out in the Full Window mode and, when 
combined together, provide a useful exposure time of 19 ks. 
We only use the EPIC pn data because of its high sensitiv- 
ity in the Fe K band. The data analysis presented here was 
performed using the latest version of the standard analysis 
package, SAS 6.0. Single and double events from the detec- 
tor were selected and the data reduction was carried out 
following the standard procedure. 
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Figure 1. The 2-10 keV band XMM-Newton EPIC pn spectrum 
of Arp220. 



3 RESULTS 

3.1 The XMM-Newton EPIC spectrum 

The hard X-ray (> 3 keV) emission in Arp220 is only re- 
solved at the resolution of the Chandra X-ray Observatory 
(Clements et al 2002; and see Section 3.2), and is point-like 
when viewed with the XMM-Newton telescope, although a 
much larger extension of soft X-ray emission is clearly re- 
solved. As we focus on the hard X-ray emission, the spectral 
data are taken from a circular region of a 25 arcsec radius, 
enough to collect most of the hard X-ray photons. The de- 
tailed location of the hard X-ray source is shown in Section 
3.2. 

The EPIC pn spectrum is shown in Fig. 1. The soft X- 
ray emission below ~ 2keV is mostly due to the extended 
nebula, within which a wide range of spectral variations are 
present between regions, as revealed by the Chandra data 
(McDowell et al 2003). The data analysis presented here is 
restricted to the energy range 2.5-10 keV to focus on the 
hard X-ray emission around the double nucleus. 

A prominent Fe Kq line is detected at 3.5a significance 
(Fig. 1). The centroid of the feature is found at 6.72±°° 5 5 
keV (in the rest frame; the errors represent the 90 per cent 
confidence region for one parameter of interest), indicating 
that FeXXV is the major line component. Fitting a single 
gaussian gives a line flux of 1.7lo!i x 10~ 6 phs -1 cm~ 2 . The 
corresponding equivalent width (EW) against the underly- 
ing continuum is 1.9±0.9 keV. The inclusion of the gaussian 
line in the fitted model reduces \ 2 by — 10 (23 degrees of 
freedom). The single gaussian fit does not require statisti- 
cally significant broadening, but a multiple line complex in 
the 6.5-7 keV range would be more likely. 

There are suggestions of other spectral features at 3.3, 
3.9, 4.2 and 5.5 keV with their detection being at 2.6a or 
lower. They all have possible identifications with a highly 
ionized plasma. Because of their low significance, we only 
give results of spectral fitting with narrow gaussians for 
those emission features in Table 1, and will not use them as 
critical materials for further discussion. However, we point 
out that they would be an important key to constrain the 
origin of the hard X-ray emission, once their detections are 
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Table 1. Emission line features in the 3—10 keV band. Results are 
obtained from fitting each spectral feature with a narrow gaus- 
sian, and the centroid energy is corrected for the galaxy redshift 
(z = 0.018). The fit with these gaussian lines gives x 2 = 1-0 f° r 
17 degrees of freedom. The errors quoted are the 90 per cent con- 
fidence region for one parameter of interest. RRC stands for ra- 
diative recombination continuum, which is only relevant for pho- 
toionized gas. tThe line flux limits are obtained by fitting with 
the energies fixed for ArXVIII and CaXX, respectively. 
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confirmed by better quality data. Especially, radiative re- 
combination continua are spectral features unique to pho- 
toionized gas (Section 4.5). 

The 2.5-10 keV continuum is very hard: fitting a 
power-law modified only by Galactic absorption (Nh = 4x 
10 20 cm~ 2 ) gives T = 1.2±° *. The total 2-10 keV flux is 
1.1 x 10~ 13 erg cm -2 s~ . This value is in agreement with 
the value obtained from the Chandra observation (Clements 
et al 2002), but is smaller than the BeppoSAX MECS value 
1.8 x 10 _11 erg cm" 2 s — l . As pointed out by Clements et 
al (2002), the MECS aperture contains two hard X-ray 
sources to the south, and the discrepancy with BeppoSAX 
is probably due to contamination from these sources. The 
corresponding 2-10 keV luminosity is 7 x 10 40 erg s _1 , and 
£2-10k C v/£ir =± 1.5 x 10" 5 . 

The possibility that these contaminating sources in the 
XMM-Newton beam emit the iron line emission can be ruled 
out. The nucleus of Arp220 is the only significant source 
within the beam in the narrow band (6-7 keV) Chandra 
image (see Section 3.2). Given the large EW of the line 
(which means the iron line emission to dominates the 6- 
7 keV band), the Arp220 nucleus is certainly the iron line 
source. 



3.2 Chandra imaging 

As reported previously, the hard X-ray emission is slightly 
extended but concentrated around the double nucleus 
(Clements et al 2002). Most of the 3-7 keV emission comes 
from radii within 0.7 kpc and very little originates beyond 
1 kpc. Fig. 2 shows the Chandra 3-7 keV band image with 
the HST NICMOS 1.6pm image superposed. The astrome- 
try of the Chandra image has been corrected using the latest 
attitude correction file, and the HST image has been regis- 
tered in the manner discussed in Scoville et al (1998). In 
this updated registration, the western nucleus (radio and 
near-infrared) now coincides with the 3-7 keV peak. Note 
the positions of the Eastern (radio) nucleus and the 1.6 fj,m 
peak, the displacement of which is presumably due to strong 
obscuration (see Scoville et al 2000 for details on the near- 
infrared extinction in this region). The centroid of the 6-7 
keV band emission, for which only 12 counts are detected 
and which is presumably mostly due to Fe K line emission 
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Figure 2. The Chandra 3-7 keV image of the double nucleus 
region of Arp220 overlaid by the HST NICMOS 1.6 ^m contour 
map (Scoville et al 2000). The positions of the double radio nuclei 
are indicated with plus symbols in green. 



(see also Clements et al 2002), may be slightly displaced 
from the Western nucleus to the East, but this is not con- 
clusive. The Chandra spectrum of the hard X-ray source is 
consistent with the spectral model with gaussian lines fit- 
ted to the XMM-Newton spectrum in the 2-5 keV range. 
At higher energies, a comparison becomes difficult as the 
efficiency of Chandra declines sharply. 



4 THE ORIGIN OF THE HARD X-RAY 
EMISSION 

4.1 X-ray binary emission 

It is generally assumed that the 2-10 keV emission in star- 
burst galaxies is dominated by integrated emission from X- 
ray binaries, and its luminosity has been found to be cor- 
related with star formation rate indicators, e.g., infrared 
luminosity (Nandra et al 2002; Ranalli, Comastri & Setti 
2003; Grimm et al 2003). There may be non-thermal inverse 
Compton scattered emission of infrared photons by relativis- 
tic electrons (e.g., Moran & Lehnert 1997), but its contribu- 
tion is probably minor. In Arp 220, the hard X-ray emission 
is resolved with Chandra and its extension coincides roughly 
with the dense molecular disk (e.g., Sakamoto et al 1999) in 
shape. However, the detection of the strong Fe K line readily 
rules out X-ray binaries as a major source of the 2-10 keV 
emission because of the spectral incompatibility. 

If the infrared luminosity from Arp220 is entirely due 
to a starburst, the implied star formation rate is w 2OOM0 
yr~ (Kennicutt 1998). With this high star formation rate, 
emission from X-ray binaries - which is estimated by as- 
suming only high-mass X-ray binaries 1 , and by following 
Franceschini et al (2003; see also Persic & Rephaeli 2002; 
Persic et al 2004) - should dominate, or in fact, exceed the 



1 In a young starburst system, as low mass X-ray binaries are 
not yet formed, it is appropriate to consider only high-mass X- 
ray binaries (Persic & Rephaeli 2002). 
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observed luminosity by more than one order of magnitude, 
even if absorption of the order of 10 22 -10 23 cm~ 2 in Nh is 
taken into account. The lack of X-ray binary emission ap- 
pears to be a general problem with a starburst interpretation 
for Arp220, given the good correlation between the 2-10 keV 
X-ray binary emission and the infrared luminosity claimed 
for star forming galaxies (e.g., Ranalli, Comastri & Setti 
2003; Grimm et al 2003; Persic et al 2004). 

In the following subsections, we discuss two possible ori- 
gins of the Fe K line in the context of a starburst, based on 
the thermal emission model for the observed X-ray spec- 
trum. Both explanations require as high as a star formation 
rate mentioned above, therefore the lack of the X-ray binary 
emission remains to be a problem. 

4.2 Thermal emission model 

With the Fe K line centred at 6.7 keV, it seems tempting 
to interpret the hard X-ray emission as thermal emission 
from hot gas associated with a starburst. Fitting the 2.5- 
10 keV data with the collisionally ionized plasma spectra 
computed by the MEKAL code (e.g., Kaastra 1992) gives 
the following results: The temperature of the gas implied 
from the fit is 7.4tg; 4 keV. The absorption column density 
is not well constrained, but the likely value for the best- 
fit temperature is 3 x 10 22 cm~ 2 . The metallicity, which is 
primarily determined by the Fe K line strength, is found 
to be 2.2lf; 2 solar (the solar abundance table by Anders & 
Grevesse 1989 is used here). 

The thermal emission model accounts for the continuum 
shape and the Fe K feature (with \ 2 = 21.5 for 17 degrees 
of freedom), but would leave the possible lower energy fea- 
tures unexplained. Since the interstellar gas in a starburst 
region is expected to be enriched by core-collapse supernovae 
(e.g., Type II SNe), which produce substantial a-elemcnts 
but a relatively small amount of iron. The twice-solar Fe 
metallicity required to explain the Fe K feature is already 
large, although it may not be surprising for a region with in- 
tense star formation (e.g., Fabbiano et al 2004). A non-solar 
abundance ratio between a-elements and Fe, as expected for 
chemical enrichment by Type II SNe, might be required, if 
the detection of the high ionization Ar and Ca features are 
confirmed by higher quality data. If confirmed, the CaXIX 
He-a at 3.9 keV (Table 1) would not be compatible with the 
single temperature thermal model, because with a temper- 
ature of 7 keV, calcium ions are mostly CaXX. 

4.3 A hot bubble in a starburst region 

According to the starburst-driven superwind model, a hot 
bubble of internally shocked wind material with a tempera- 
ture of several keV will form in the starburst region (Cheva- 
lier & Clegg 1985), which eventually breaks away to drive a 
galactic-scale outflow (e.g., Tomisaka & Ikeuchi 1988). Al- 
though such a hot bubble is expected not to be radiative be- 
cause of its rarefied interior, and is indeed rarely observed di- 
rectly in starburst galaxies (e.g., Suchkov et al 1994; Strick- 
land & Stevens 2000; Hoopes et al 2003) 2 , the temperature 

2 The presence of such hot gas has been reported in M82 by 
Griffiths et al (2000) and in NGC253 by Pictsch ct al (2001), but 



of kT « 7 keV implied from the above thermal emission 
model matches the prediction. In fact, the observed temper- 
ature and luminosity agree roughly with the predicted values 
from the wind solution by Chevalier & Clegg (1985) with the 
mass and energy input rates estimated by Starburst99 (Lei- 
therer et al 1999) for the star formation rate of 200 Mq yr -1 
(see below) and a high thermalization efficiency (~ 1). 

As the Chandra image shows, most of the hard X-ray 
emission in Arp220 comes from within « 0.7 kpc of the 
double nucleus. Assuming the thermal plasma is uniformly 
distributed within a sphere with a radius of 0.7 kpc, the 
mean gas density and the total gas mass are derived to be 
n gas ~ 0.6 cm~ 3 and M gas ~ 2 x 10 7 M Q , respectively. The 
thermal energy contained in the volume is estimated to be 
~ 3.5 x 10 erg. The bolometric luminosity of the gas of 
« 1.6 x 10 41 erg s _1 means that the cooling time of the gas is 
~ 7 x 10 7 yr. This value is insensitive to the non-solar abun- 
dance ratio, because, at a temperature of kT ~ 7 keV, the 
cooling is dominated by bremsstrahlung continuum rather 
than line emission. 

With the star formation rate of 200 Mq yr -1 , the mass 
injection rate from OB stars is ~ 3OM yr -1 , assuming 
a Salpeter IMF and a mass loss rate of 10 -5 Mq, typical 
for each massive star (see also Elson, Fall & Freeman 1989; 
Heckman et al 1990). Starburst99 also gives a mass injection 
rate of ~ 5OM yr -1 both from stellar winds and super- 
novae, assuming continuous star formation and a starburst 
age > 10 7 yr. The above IMF is assumed to have an upper 
mass cut-off at 100 Mq. If the IMF is truncated at 30 Mq, 
the estimate is reduced by a factor of 2 or less. However, 
given the cooling time, the expected mass input is sufficient 
to supply and maintain the hot gas. 

At a temperature of 7 keV, the sound speed is c s ~ 800 
km s~ . So, the sound-crossing time over the radius of the 
hard X-ray emitting region (0.7 kpc), is 8 x 10 5 yr. Since 
the sound crossing time is much shorter than the cooling 
time, the hot gas, which will not be static because of its 
high pressure (~ 10~ s dyne cm -2 ), will drive a superwind. 

In terms of the energetics of the gas, heating by SNe is 
sufficient to counterbalance the radiative loss. As discussed 
in the following section, the expected supernova rate is ~ 2 
SNe yr -1 , implying an energy injection rate of ~ 6 x 10 43 
erg s _1 (or Starburst99 gives the mechanical luminosity of 
~lx 10 44 erg s _1 ). The required fraction of the SN energy 
that goes to heating is very small rjh ~ 3 x 10~ 3 . Therefore, 
the bulk of the energy has to escape in the form of a wind - 
given the moderate soft X-ray nebula luminosity - without 
depositing the energy onto the galactic medium to radiate. 
This hot bubble interpretation appears entirely plausible, 
apart from the X-ray binary problem mentioned above. 

4.4 Luminous radio supernovae 

An alternative origin for the thermal emission is an ensem- 
ble of supernovae. Several radio knots distributed over the 
nuclear region have been imaged with VLBI, and are inter- 
preted as very luminous compact radio supernovae taking 
place in a dense environment (Smith et al 1998). They can 

for NGC253, photoionization by AGN is proposed by Weaver et 
al (2002). 
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be luminous X-ray sources, as Type Iln SNe (Benetti et al 
1995) like SN1986J (Houck et al 1998), SN1988Z (Fabian & 
Terlevich 1996) and SN1995J (Fox et al 2000) have X-ray 
luminosities of 10 40 -10 41 erg s _1 . Their X-ray spectra seem 
hard enough to match the hard-band spectrum of Arp220, 
i.e., temperatures of 3-10 keV in kT when fitted with a 
thermal emission model (as summarised in Fox et al 2000), 
and at least SN1986J shows evidence for a strong Fe K line 
at 6.7 keV (Houck et al 1998). With a star formation rate 
of 200 Mq yr _1 , the expected supernova rate is ~ 2 SNe 
yr~\ as also estimated by Smith et al (1998). If 1 per cent 
of the total energy of each SN goes into radiation, and is 
emitted at 10 41 erg s" 1 , then the cooling time is ~ 3 yr, 
which is roughly comparable to the estimated average time 
between SNe (however, as the radiative efficiency might be 
much higher in the dense environment, these SNe could over- 
produce the integrated X-ray luminosity). Therefore it is 
plausible to have multiple X-ray SNe at the same time, and 
to maintain a stable hard X-ray luminosity until most of the 
massive stars die out. 

Whether powerful SNe like SN1986J are constantly pro- 
duced in the Arp220 nuclear region is questioned by the 
second VLBI observation following the results presented in 
Smith et al (1998) three years later. Most of the radio knots 
have faded by 8-50 per cent (some remain at the same 
brightness), but no new compact sources have appeared 
(Smith et al 1999). The slow decline is not consistent with 
a luminous SN model. The SN rate also has to be lowered 
(< 0.3 SNe yr" 1 , Smith et al 1999). Although this esti- 
mate is only for the powerful radio SNe, X-ray luminous SNe 
are probably closely related. With these facts, it is unclear 
whether the above explanation of the hard X-ray emisson 
with powerful SNe is sustainable. 

4.5 Photoionized gas 

Photoionized gas illuminated by a hidden AGN is still a vi- 
able alternative for Arp220. In fact, if the possible radiative 
recombination continuum (RRC) feature is real, it may be 
the most likely interpretation. A high ionization parame- 
ter i = L lon /{nR 2 ) ~ 10 3 erg cms 1 , where Li an erg s 1 , n 
cm~ 3 , R cm, are respectively the ionizing luminosity, density 
and distance from the source, is implied by the line emission 
centred at 6.7 keV. The extension of the hard X-ray emis- 
sion points to the presence of a larger region of low density 
interstellar gas with density of the order of 1-10 cm -3 . Ex- 
tended photoionized nebulae have been found in a number 
of nearby Seyfert 2 galaxies, some of which show extended 
Fe K emission (e.g., NGC4945, Done et al 2003; NGC4388, 
Iwasawa et al 2003). 

The photoionized gas in AGN appears to have ioniza- 
tion parameters distributed over a wide range (e.g., Krolik 
& Kriss 2001; Kinkhabwala et al 2002). Although detailed 
modelling of the Arp220 hard X-ray spectrum is beyond the 
scope of this paper, given the quality (and also the spec- 
tral resolution) of the present data, our preliminary inspec- 
tion using the photoionization code, XSTAR (Kallman & 
Bautista 2001) indicates log £ in the range of 2.8-3.5 is rel- 
evant to the observed spectral features. Absorption with a 
column density of the order of 10 23 cm~ 2 would be required 
for this photoionized spectrum not to dominate the soft X- 
ray band and to make the spectrum as hard as is observed. 



This photoionization model predicts RRC features 
at 5-5.5 keV from CaXIX,XX and at ~ 9 keV from 
FeXXV,XXVI, which should be relatively isolated so that 
they could be resolved at the CCD resolution if the data 
quality is improved. This could provide a crucial test for 
the photoionization model when longer exposure hard X-ray 
data are obtained. 

Strong Fe K line complex emission, consisting of a 6.4 
keV line, which is usually the strongest as expected from 
Fel-XVII, and higher energy lines of FeXXV (6.7 keV) and 
FeXXVI (6.97 keV), has been observed in a number of 
Compton-thick Seyfert 2 nuclei (e.g., NGC1068, NGC6240, 
NGC4945). It is clear that the 6.4 keV line and the higher 
energy lines originate from different matter because of their 
large difference in ionization stages. The narrow line width 
of the latter naturally suggests that the line emitting gas is 
optically thin. 

The strong 6.4 keV lines seen in good quality spectra 
of Seyfert 2 nuclei often show a Compton shoulder, i.e. a 
weak redward extension to the line core (e.g., Iwasawa et al 
1997; Bianchi et al 2002; Kaspi et al 2002; see also George 
& Fabian 1991; Matt 2002 for theory). This means that the 
6.4 keV line results from reflection from optically thick mat- 
ter, which is, in the torus model for the unification scheme, 
identified as the visible surface of the inner wall of the ob- 
scuring torus (e.g., Awaki et al 1991; Ghisellini, Haardt & 
Matt 1994; Krolik, Madau & Zycki 1994). 

No detection of a 6.4 keV line in the Arp220 spectrum 
implies the lack of reflection from optically-thick cold mat- 
ter, and perhaps suggests that the visibility of the torus in- 
ner wall is very limited. The conditions for an energetically 
significant AGN to exist in Arp220, imposed by the hard X- 
ray observation, are that the obscuring matter must have a 
column density larger than Ah = 10 25 cm~ 2 with a covering 
factor close to unity (Iwasawa et al 2001). With the very 
high covering factor, the inner wall of the obscuring torus 
is capped and never visible, which could explain the lack of 
a 6.4 keV line. If a tiny fraction of light from a hidden nu- 
cleus leaks through the heavy obscuration, the surrounding 
low density medium could be highly ionized to give rise to 
"hot" Fe K line emission. 

The covering fraction, (1 — /), can be estimated based on 
the Fe K line flux, by comparing the well-studied Compton- 
thick AGN, NGC1068. NGC1068 has a large infrared excess 
with Lir ~ 7 x 10 44 erg s _1 for an assumed source distance 
of 14.4 Mpc. The 6.8 keV line flux measured in the Chandra 
data is 5.5 x 10 _5 phs _1 cm~ 2 (Young, Wilson & Shopbell 
2001; a similar flux was measured by Iwasawa et al 1997 
for the sum of FeXXV and FeXXVI in the ASCA data). 
Kinkhabwala et al (2002) estimated fLi 0n ~ 10 43 erg s _1 , 
by fitting the detailed photoionization model to the RGS 
data on NGC1068. If both NGC1068 and Arp220 are pow- 
ered by hidden AGN, and the hot Fe K line emitting region 
in NGC1068 sees the same ionizing source, then fLi 0n for 
Arp220 is ~ 0.8 x 10 43 erg s"\ and with L lon ~ (l/2)L ir , 
/ is estimated to be ~ 0.4 per cent, which can be com- 
pared with / ~ 3 per cent for NGC1068 (/ < 10 per cent 
from the Chandra image, Kinkhabwala et al 2002) under 
the same assumption. Thus, reducing the opening fraction 
of the obscuration further would make the photoionization 
model consistent with that for NGC1068. 
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5 THE POWER SOURCE OF THE 
FAR-INFRARED EMISSION 

While there is no doubt that a starburst is taking place 
in Arp220, various mid-infrared characteristics, summarised 
by Spoon et al (2004, and reference therein), the superwind 
luminosity inferred from the Ha and soft X-ray nebulae 
(Suchkov et al 1996; Iwasawa 1999; also McDowell et al 2003 
for an alternative interpretation of merger shock), and the 
X-ray binary emission in the hard X-ray band, as discussed 
above, are all unusually low relative to the far-infrared lu- 
minosity for starburst galaxies. 

The peculiarity of Arp220 is well illustrated by Spoon 
et al (2004), who re-examined the ISO mid-infrared spec- 
trum with the knowledge of a number of ice absorption 
features across the band. Their spectral analysis and the 
higher resolution study by Soifer et al (2002) demonstrate 
that the mid-infrared emission in Arp220 consists of diffuse 
PAH emission with a moderate amount of absorption and 
a NGC4418-like, heavily absorbed continuum. The latter is 
more likely to contribute to the far-infrared luminosity of 
Arp220, although the above authors note that, if all the far- 
infrared luminosity is powered by a hidden AGN, it would 
not even be visible at mid-infrared wavelengths. 

If, say, 10 per cent, of the bolometric luminosity (~ 
lO 11 ^©) was due to a moderately obscured starburst, the 
outward characteristics of such a starburst as described 
above would then be in agreement with most other starburst 
galaxies. For example, the soft X-ray nebula has the lumi- 
nosity of ~ 1 x 10 41 erg s _1 , consistent with the predicted lu- 
minosity ratio of X-ray and starburst-powered infrared emis- 
sion for a superwind nebula (log (SX/IR) « —3.5, Leitherer 

6 Heckman 1999; Strickland & Stevens 2000). Also, the X- 
ray binary emission in the hard X-ray band approaches the 
correlation line with the star formation rate for local star- 
burst galaxies (Ranalli et al 2002; Grimm et al 2002). 

If such a lower luminosity starburst was indeed the case, 
the rest of Arp220's luminosity would have to be explained 
by something else. Spoon et al (2004) suggested a deeply 
shrouded ultra dense starburst, in addition to the less ob- 
scured starburst, which is responsible for the diffuse PAH 
emission. The estimated optical depth based on their mod- 
elling of the ISO spectrum implies that the column density 
to such a star cluster would be AT H ~ 10 23 cm~ 2 . Since X- 
rays at energies higher than a few keV are transparent to 
this level of obscuration, luminous hot gas, for instance, as- 
sociated with the star cluster would be observable as a hard 
X-ray excess, which is not apparent in the data. 

An alternative option is an even more deeply buried 
AGN. If the photoionized gas interpretation was correct for 
the detected Fe K line, then the presence of an AGN would 
be proved. It is not straightforward to estimate the intrin- 
sic luminosity of such an obscured AGN without seeing its 
transmitted radiation: once the line of sight Thomson depth 
exceeds unity (Compton thick), the visibility of nuclear X- 
ray emission is largely determined by geometry, e.g., the cov- 
ering factor, distribution of the gas density etc., which are 
not known. However, the comparison with NGC1068 (Sec- 
tion 4.4) demonstrates that with / ~ 0.4 per cent, the de- 
tected Fe K line luminosity is consistent with the assumption 
that all the luminosity is powered by an AGN. Therefore 
once the photoionization origin for the Fe K line is con- 



firmed, it would be plausible that an AGN dominates the 
energetics. The radio spectrum of Arp220 must then be ex- 
plained by free-free absorption, the opacity of which should 
be larger than those in the Compton-thick AGN in NGC6240 
and NGC4945 (see Fig. 5 in Iwasawa et al 2001). 

There are a few interesting infrared luminous objects 
to compare with Arp220. NGC4418 is a luminous infrared 
galaxy with L\ T ~ 1 x lO 11 ^©, and its mid-infrared spectrum 
shows no PAH features, but is shaped by silicate and various 
ice absorption features (Spoon et al 2001). Very faint X-ray 
emission with Lx ~ 10 39 erg s _1 has been detected from the 
centre of NGC4418 with Chandra (Maiolino et al 2003). The 
luminosity ratio of X-ray and infrared emission for NGC4418 
is even smaller than that of Arp220. The hyper-luminous in- 
frared galaxy, IRAS 00182-7112 (z = 0.327), has kinemat- 
ical evidence for a superwind (Heckman et al 1990), as in 
Arp220, and its mid-infrared spectrum is classified as Class 
1 by Spoon et al (2002), i.e., it is NGC4418-like. A recent 
XMM-Newton observation of IRAS 00182-7112 detected a 
Fe K line at 6.7 keV with EW ~ 1.5 keV on a flat continuum 
(K. Nandra, priv comm), reminiscent of Arp220. 

These three objects share similar properties, and if they 
are powered by stars, they are starbursts for which the 
widely used relation between far-infrared (or the translated 
star formation rate) and hard X-ray (X-ray binary emis- 
sion) luminosities does not apply. Their 2-10 keV X-ray lu- 
minosities are well below the correlation line of Ranalli et al 
(2003) and Grimm et al (2003). The detection of the Fe K 
line, in Arp220 in particular, means that their X-ray binary 
contribution to the 2-10 keV band must be minor, mak- 
ing them deviate from the correlation further. Hiding most 
of the X-ray binaries behind thick obscuration required a 
substantial mass of gas and dust, given the observed spa- 
tial extension (e.g., Arp220). If the formation of X-ray bina- 
ries were suppressed, an unusually large suppressing factor 
(~ 20 for Arp220 when compared with Grimm et al 2003) 
would be required. An heavily obscured AGN may offer a 
more relaxed solution for the embedded power sources for 
all three of these objects. For Arp220, of course, a substan- 
tial starburst (~ 10 11 Lq) would still be required to explain 
the near- mid infrared light, superwind signatures and radio 
supernovae etc. 
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